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Abstract 

k>( ' It has recently been shown that the maximal kinematical invariance 

group of polytropic fluids, for smooth subsonic flows, is the semidirect 
product of SL(2,R) and the static Galilei group G. This result purports 
to offer a theoretical explanation for an intriguing similarity, that was re- 
cently observed, between a supernova explosion and a plasma implosion. 
In this paper we extend this result to discuss the symmetries of discontin- 
uous flows, which further validates the explanation by taking into account 
shock waves, which are the driving force behind both the explosion and 
implosion. This is accomplished by constructing a new set of Rankine- 
Hugoniot conditions, which follow from Noether's conservation laws. The 
new set is dual to the standard Rankine-Hugoniot conditions and is related 
to them through the SL(2, R) transformations. The entropy condition, 
that the shock needs to satisfy for physical reasons, is also seen to remain 
invariant under the transformations. 
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1 Introduction 



It has recently been observed that the density profiles of a supernova explosion 
and an inertial confinement plasma implosion^ [2CU are strikingly similar. An 
empirical basis for this intriguing duality between explosion and implosion was 
given by Drury and Mendonga |2J who pointed out that Euler's equations of fluid 
dynamics, which describe both the systems, are form-invariant under a set of 
nonlinear coordinate transformations viz. x — > x/t, t — ► — 1/t. The minus sign 
in the time transformation maps an explosion to an implosion and the inversion 
allows large time scales to be mapped to small time scales and vice versa. These 
transformations suggest that the maximal kinematical invariance group Q of 
fluid dynamics is larger than the standard Galilei group. It is now known that 
this larger group is a twelve-parameter semidirect product, Q = SL{2, R) A G 
jSJ|H], where G is the nine-parameter, connected, static Galilei group: 

x — > Rf + vt + a, t —> t (1) 

and SL(2, R) is the group consisting of the transformations: 

cti.4- ft r 
t -» —4 , x -» — — with a6-py = l. (2) 

Physically, the three-parameter SL(2, R) group consists of time translations, 
scale transformations, and a one-parameter set of time-dependent scale trans- 
formations called expansions. The transformations proposed by Drury and Men- 
donga are a special case of the SL(2, R) transformations with (a, (3, 7, S) — 
(0, — 1, 1, 0). The SL(2,R) part of Q is therefore important for a better 
understanding of the explosion-implosion map. 

It should be pointed out that the naive expectation of using time-reversal in- 
variance, to explain the similarity between explosion and implosion, is untenable 
here since the length and time scales involved in the two systems are drastically 
different. Invoking scaling arguments is not of much help since, although a 
composition of time-reversal and suitable scalings leaves the equations of fluid 
mechanics and the Reynolds number invariant, it has the property of reversing 
the direction of time's arrow and thereby violates the second law of thermody- 
namics. As a consequence, when applied to a shock wave, such transformations 
violate entropy conditions that define the physicality of the shock. As is well- 
known, however, both the supernova explosion and the plasma implosion are 
driven by the formation and propagation of a shock wave. It is therefore impor- 
tant to examine whether the SL(2, R) symmetry, that purports to explain the 
observed duality, respects the physicality of the shock wave. With this in view, 
we extend the study of (SJ - in which the explanation of explosion-implosion du- 
ality based on the symmetry group Q was restricted to smooth, subsonic flows 
- to examine shock waves. It will be shown that shock wave solutions are con- 
sistent with the symmetries of the maximal kinematical invariance group Q in 
the following precise sense. 

A shock in a fluid is described mathematically by the well-known Rankine- 
Hugoniot jump conditions [Jj. So the natural question to ask is: What happens 
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to these conditions under the action of the SL(2, R) group? This question is 
best answered not in the framework of the partial differential equations of fluid 
dynamics, but by reverting back to their so-called primitive form i.e. expressing 
them as conservation laws. The conservation laws are completely equivalent to 
the partial differential equations for smooth flows, but produce the Rankine- 
Hugoniot jump conditions for discontinuous flows in a natural and well-defined 
manner. The connection with SL(2,R) is made by appealing to Noether's 
theorem which asserts that corresponding to every continuous symmetry, there 
exists a conserved charge. 

Anticipating that the SL(2,R) transformations will mix the conservation 
laws corresponding to various symmetries, we construct the Noether charges 
corresponding to them and the boost transformations, in addition to the well- 
known ones for rotations, space and time translations. We then use the atten- 
dant conservation laws to establish a new set of jump conditions. It turns out 
that the new conditions are identically satisfied if the standard conditions for 
mass, momentum and energy conservation are satisfied. Although seemingly 
redundant because of this reason, the new set holds independently; following, 
as it does, from the symmetries of the fluid equations. In fact, these conditions 
are useful to prove the form-invariance of the Rankinc-Hugoniot conditions un- 
der the SL(2, R) transformations. Thus, to each physical system governed by 
the fluid dynamics equations two independent, but physically equivalent, sets 
of jump conditions can be associated, the two being related by SL(2, R) trans- 
formations. We conclude that the SL(2, R) transformations map the Rankine- 
Hugoniot conditions of the explosion to the dual Rankinc-Hugoniot conditions 
of the implosion and vice versa. Further, by specialising to the Drury-Mendonga 
transformations, x — > x/t, t — > —1/t, we show that the jump conditions for 
boosts and expansions, along with the continuity equation for mass conservation, 
provide an independent, albeit equivalent, description of the shock. They may 
be viewed either as the dual of the standard Rankinc-Hugoniot conditions, or, 
in the language of passive coordinate transformations, as the standard Rankinc- 
Hugoniot conditions in the dual coordinate system corresponding to the choice 
(a, /3, 7,5) = (0,-1,1,0). Similar dual conditions exist for each choice of the 
SL(2, R) parameters. 

It is well-known that Rankinc-Hugoniot conditions describe not only shocks, 
but other discontinuities like slip and contact discontinuities, for example. There- 
fore, the map between the dual sets of Rankinc-Hugoniot conditions would be 
relevant to explosion-implosion duality only if both the sets refer to shocks. In 
other words, only those Rankinc-Hugoniot conditions that describe a shock and 
only those SL(2, R) transformations which map a shock to a shock are of in- 
terest for explosion-implosion duality. Moreover, Rankinc-Hugoniot conditions 
say nothing about the physicality of the shock - this information is contained 
in additional inequalities for its entropy that the shock needs to satisfy. Physi- 
cal shocks are distinguished from others because their entropy always increases 
across the shock front. We verify explicitly that this requirement is unaffected 
by the SL(2, R) transformations. 

The physicality of the map between explosion and implosion may also be es- 
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tablished in the following subtle manner: Although the notions of viscosity and 
heat conduction lose their meaning in the immediate vicinity of the shock, be- 
cause the changes in all the quantities they depend on are so great, they do play 
an important role in the formation and maintenance of a shock discontinuity 
In particular, the positivity of the coefficients of viscosity and heat conduc- 
tion guarantees that the shock satisfies the appropriate entropy conditions El- 
Hence, Euler's equations ought to be considered as a special case of a more 
general set of fluid equations with vanishingly small viscosity. Requiring the 
sign of the viscosity to remain unchanged under the transformations establishes 
the physicality. The Navier-Stokes equations - which are the obvious choice for 
including viscosity - are not invariant under the full SL(2, R) part of Q, but 
only under the standard Galilean transformations. However, a more general 
set of fluid equations with viscosity fields transforming appropriately under the 
SL(2, R) transformations has a maximal kincmatical invariance group given by 
Q |H] • Hence we use these equations for our purpose of examining the behaviour 
of non-vanishing viscosity under the SL(2,R) transformations. Similar argu- 
ments apply for heat conduction, but it does not bring in any new qualitative 
features and hence is omitted from further discussion. 



2 Symmetries of Fluid Dynamics 

In this section we briefly recapitulate the results of [5]. The general fluid equa- 
tions in n-dimensional space are [0] 

Dp = -pV-u (3) 
pDu = -Vp + V (4) 
De = -{e + p)V-u (5) 

where 

d 

D = -+u-V 
at 

and 
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Vi = Vj yi(\7jUi + ViUj - -5 i3 V k u k ) J + Vi(CVfcUfc) 

In the above equations p, u, p, e stand for the density, velocity, pressure and 
energy density of the fluid respectively and ?y, £ are the bulk and shear viscosity 
fields. 

The above differential equations are usually augmented by an algebraic con- 
dition called the polytropic equation of state which relates the pressure to the 
energy density as 

P = (7o - l)e (6) 

where 70 is a constant called the polytropic exponent. As shown in [3], the 
maximal invariance group of the above set of equations is Q = SL(2, R) A G, 
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provided the polytropic exponent takes the standard value for an ideal, nonrel- 
ativistic fluid viz. 70 = 1 + — . For this value, the fluid equations are invariant 
under the following transformations [S]: 

Connected, static Galilei transformations: 

Let g denote a general element of this sub-group then 

g : t' — t , x = Ric + vt + a (7) 

with R an orthogonal matrix. Under the action of g, the fields p and u transform 
as 

p 1 = p and u' = u + v (8) 
SL(2,R) transformations: 

Let cr denote a general element of the SL(2, R) part of Q then 

, at + P _. x 

a : t = , x = where ad — #7 = 1 (9) 

jt + 6 jt + S 

Under the action of cr, the fields transform as 

p = (jt + S) n p and u' = (7* + S)u - jx (10) 

For both g and cr, the transformations of e and p can be worked out once the 
transformation properties of p are known since 

e = xp 7o i p = (7o _ 1)e t (n) 

with the field x ~ related to entropy - transforming like a scalar. The transfor- 
mation properties of the viscosity fields are similar to the density p 

rf = (lt + 6) n ri and C'=(7* + <5)"C (12) 

The above results were derived in jS] by requiring the invariance of the Action 
for the simple case of an inviscid and isentropic fluid. The symmetry of the 
equations followed by subsequently relaxing the simplifications to arrive at the 
general fluid equations, ft should be noted that the requirement of the invariance 
of the Action is sufficient, but not necessary, for the form invariance of the equa- 
tions that follow from it. Any transformation that leaves the Action invariant 
upto a multiplicative factor produces equations of motion which have the same 
form. If this is taken into account, the condition aS—fij = 1 is no longer required 
and SL(2, R) gets replaced by GL(2, R) in the maximal invariance group of the 
general fluid equations. However, it is sufficient for our purposes to concen- 
trate on the variational symmetries of the fluid equations and for this purpose, 
Q = SL(2, R)AG. Finally, for the sake of completeness, it should also be pointed 
out that the SL(2, R) condition is invariant under the following discrete symme- 
tries (a, (3, 7, 5) -» (a, -/?, -7,5), -> (~a, (3, 7, -6), -» (-a, -f3, -7, -S) 
of the SL(2, R) parameters. 
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3 Conservation Laws 



In this section we construct the conservation laws corresponding to the symme- 
tries outlined in the previous section. In order to do this, it is useful to revert 
back to the Action formalism and obtain the results for the subclass of inviscid, 
isentropic and irrotational flows. The corresponding expressions for a general 
fluid can then be worked out along the lines of |Sj ■ 

For inviscid, isentropic and irrotational flows, the Lagrangian density is given 

by 

C = p{^-\{V<}>f)~p<° (13) 

where V</> stands for the curl-free part of the velocity vector field u. Let 
/i = 0,1,2,3 and x^ be a four- vector under the transformations of the pre- 
vious section i.e. x l with i — 1, 2, 3 are the components of x and x° — t. Let 
the infinitesimal variations in the coordinates and fields be defined as 

5x» = x^' - x" and 5<j>{x) = <j>'(x') - cf>(x) (14) 

Then the variations for translations, rotations, boosts, dilatations, and expan- 
sions respectively are given by 

5x t * = a tl , 8x l = x ] , 6x l = v l t, 5x i = Xx\ St = 2Xt and Si* = -ptx^ 

(15) 

where the parameters A, a , [i are expressible in terms of the SL{2, R) param- 
eters a, (3, 7, 5. The field variation is given by 

The variation in p is not important since no derivatives of p appear in the 
Lagrangian density. Using these results we find, by a straightforward application 
of Noether's theorem jTH], that the following quantities, integrated over all space, 
are constants of motion: 

Temporal Translations : H = ^(V0) 2 + p la (17) 

Spatial Translations : P = pVcj) (18) 

Rotations : L = P x x (19) 

Boosts : K = Pt- px (20) 

Dilatations: D = -2tH + x ■ P (21) 

Expansions : A = t 2 H - tx ■ P + ^x 2 (22) 
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The conditions of irrotationality and isentropicity can be relaxed easily and one 
sees that Euler's equations 

p = -V • (pit) (23) 



pit= -p(u-V)u-Vp (24) 

e = -V-(£u)-pV-u (25) 
can be expressed in the form of conservation laws, 

d d , 

at^-a^ (26) 

— (ptii) = - — {puiUj + Sijp) (27) 

|(^ 2 + e) = -A[(l^ +e+pK] (28 ) 

for mass and the translation generators found above. These can be reexpressed 
succinctly as follows: 

Mri = ' M% = 0, and M% = (29) 

The zeroeth components of the above currents, namely p, pu, and ^pw 2 + e, give 
the charge densities which, when integrated over all space, give the conserved 
charges. As is well-known, these are merely statements of mass, momentum 
flux, and total energy conservation. The corresponding current densities are 

4 = puj (30) 



Jp t = pu^j + S^p (31) 

J j H = { l -pu 2 +E+p)u j (32) 

The conservation laws corresponding to rotations, boosts, dilatations and ex- 
pansions can be stated similarly 

d^ £) =0, 9^=0, d^Jfa, and d, Jfc = (33) 

The charge densities are shown in (19) -(22) respectively, and the corresponding 
currents are 

Jhi = tikiXkJpt (34) 

JKi = tJp z - XiJp (35) 
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J o = Xi Jp i — 2tJn (36) 

J a = \x 2 J p - tx.J P% + t 2 J H (37) 

It may be mentioned that the above results are not surprising in the light of 
where corresponding results for a free, nonrelativistic, point particle were 
obtained through a discussion that essentially parallels the above. The note- 
worthy linear relations between the currents will, however, play a crucial role in 
this paper when we consider flows with discontinuities. 



4 Discontinuous Flows and Jump Conditions 

As long as the flows are smooth, i.e. the functions p,u,p,s G C 1 in their 
dependence on x and t, the systems (23 - 25) and (26 - 28) are equivalent. 
However, real flows are not always smooth and can develop discontinuities as 
time elapses. Such flows are described by weak solutions of differential equations 
[7] . A weak solution is generally piecewise smooth. The smooth parts satisfy the 
differential equation in the usual, or strong, form, but that does not generally 
suffice to determine the course of motion for initial data, and the equation must 
be supplemented by jump conditions. The resulting jump conditions are most 
clearly derived from the conservation laws. 

By definition any, possibly non-smooth, function J^(x, t) that satisfies 

J d^wix, t) J»(x, t) d 3 x dt = (38) 

for all test functions w(x,t) is said to be a weak solution of the differential 
equation d^J^ = 0. 

We now use the above definition to obtain the jump conditions associated 
with the system of conservation laws derived in the last section. Suppose J M (a;, t) 
has a jump discontinuity across a hyper-surface S in x, t space, while otherwise 
being continuously differentiable in some neighbourhood J\f of S (see Fig.l). Let 
w(x, t) be a test function with support in J\f. Let 1Z be the part of the support 
of w(x,t) that lies on one side of 5, say the right. Then, by Gauss's theorem 

/ df,wJ^d 3 xdt+ I wd^d 3 xdt= I df,(wjv) d 3 xdt = [ wn^J* d 3 S 
Jtz Jtz Jk Js 

(39) 

since w(x, t) = on the boundary of 1Z except on S. The second integral in the 
above equation is zero, because the conservation law holds in the strong sense in 
the interior of 1Z. Here, n(x, t) is the outward normal vector to the hypersurface 
iS. Therefore, if we integrate similarly over the left part of the support of w, 
add the result and make use of JSSJl, we find that: 

= / wn^AJ^dS (40) 
Js 
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Figure 1: Diagram for the jump condition. 



where A/ denotes the difference of the two limiting values of a function / on 
the two sides of the hypersurface S i.e. the jump of the function. This result 
follows because the vector n^, which by convention points outwards, flips its 
sign on the left side of the support. Since w is an arbitrary test function, the 
above equation implies the jump condition 



AJ M = on S . (41) 



Applying (|41|) to the conservation laws (29) for , J^p and J(jj) , we obtain 

= n^AJfo , (42) 
= rvAJgjj , (43) 
= n„AJf ff) (44) 

From here the standard Rankine-Hugoniot conditions can be derived in their 
usual form [7|. Similarly, one can apply (|41|l to the conservation laws (33) for 
Jf_ , , Jr_., and J^.s to obtain a new set of jump conditions: 

(L) (K) (D)' (A) r 

= r^AJ^ (45) 
= n.AJ^ (46) 
= ?V A J,% (47) 



? A 7 M 



= n„Aj£.y (48) 



Since the coordinates x and i are continuous on S, these conditions are all 
identically satisfied because of the jump conditions for mass, momentum and 
energy conservation, in (42-44) - a fact that can be easily verified using (19-22) 
and (34-37). 
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5 The Dual Rankine-Hugoniot Conditions 



We have seen that the new set of jump conditions associated with rotations, 
boosts, dilatations, and expansions, follow from the jump conditions associated 
with mass, momentum and energy. This suggests that the Rankine-Hugoniot 
conditions are invariant under the full kincmatical invariance group of smooth 
flows, including the SL(2,R) part. 

To see this explicitly, we consider the transformation properties of the con- 
served currents under SL(2,R). Let us begin by considering the simplest of 
these, namely the time-component of J^ p y i- e - P- From equation (10), now with 
n = 3, 

p' = (t* + 5fp (49) 

The (■yt + S) 3 factor is cancelled by an identical factor coming from the change 
of variables when we perform an integration over all space. Moreover, the trans- 
formation does not mix p with any other current. Thus, p transforms under the 
singlet representation of SL(2,R) as a scalar density. Let us now consider the 
transformation of the time-component of J^py i- e - P = pu. From (10) it now 
follows, after a little algebra, that 

P' = p'u' = {-ft + S) 3 (SP + -fK) (50) 

Thus the transformation of the spatial translation generator mixes it with 
the boost generator together with which, it forms a doublet representation of 
SL(2, R), with the prefactor (jt + S) 3 now making it a vector density. The latter 
fact is, in fact, generic to the time-components of all the currents. Likewise, we 
may consider the generator of time translations, namely the Hamiltonian, and 
it follows that 



H' = (-ft + S) 3 (j 2 A- S^D + S 2 H) (51) 

Thus the transformation of the time translation generator mixes it with the 
generator of dilatations and expansions, the three of them form the triplet (or 
adjoint) representation. The transformation properties of the rest of the cur- 
rents can be similarly worked out and the results summarised as follows: If the 
(abstract) symmetry generators T r transform as 

T; = o^TrO = M rs (a)T s , (52) 

s 

where the matrix M(a) is determined by the group structure of SL(2, R) A G, 
then the corresponding currents transform as 

/ Bt \ Bt^' « , 

W) =det b h^E M rs (a)J:(x) . (53) 
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Assembling the currents in a column, 



J (P) 
V J (V 



(54) 



one has for the transformation matrix, 



M 



/l 














\ 





a 
















7 


6 




















a 


-a(3 


/3 2 











—2aj 


(/?7 + a<5) 


-2/3(5 


\o 








7 2 


— 7<5 





(55) 



Using a5 — P"f — 1, and the fact that the determinant of a block diagonal 
matrix is the product of the determinants of the blocks, it is easy to check that 
the matrix M has unit determinant. As already pointed out, the fact that the 
currents transform like vector densities is reflected in the temporal components 
picking up a multiplicative factor (jt + <5) 3 . The spatial components follow the 
GXcimpiG 

J( P) = (it + S) n+1 r {p) - l^ht + S) n J ( ° p) (56) 

with the same 5X(2, R) transformations defined by the matrix M. 

The dual Rankinc-Hugoniot conditions are now easily obtained. The normal 
vector appearing in the jump condition (|41|l transforms like a covector, 



(57) 



so the transformed jump condition for J r is 

dx 



n'AJt?' oc det 



^2 Mrs{o)n M AJ?(x) = on S . 



(58) 



Since the determinant is smooth across the surface S, the factor in front of the 
sum can be omitted. The transformed jump condition is therefore a linear com- 
bination of the original jump conditions. In particular, the conditions for J( p ), 
Jtp\ and J(h) (the Rankine-Hugoniot conditions) become linear combinations 
of the jump conditions for J( P ), J(py ^(ity J(d) an d J(A)> 



K AJ ( P ) K " 



n' AJ^i, oc n 



(P 



A / p 

,(5 2 AJ» H) 



7 2 AJ ( %). 



(59) 
(60) 
(61) 
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The standard Rankine-Hugoniot conditions (42-44), in conjunction with the new 
set of jump conditions (45-48), then imply that the right hand side of the above 
equations is identically zero i.e. the Rankine-Hugoniot conditions are form- 
invariant. In particular, this holds for the Drury-Mendonga transformation 
t — > —1/t, x — > x/t used to relate the explosion and implosion problems. For 



this, (a, fi, 7, S) — (0, —1, 1, 0) and it follows that 

= n M Aj£) , (62) 

= nil A(xi Jfc-tJfa), (63) 

= rvA(-t 2 Jfa + taiJfa - \x 2 Jl } ) . (64) 



where we have substituted the explicit expressions for the currents J'V and 
J(Ay The conditions (62-24) are the dual Rankine-Hugoniot conditions. If an 
explosion is described by the standard Rankine-Hugoniot conditions, the corre- 
sponding implosion, obtained by a Drury-Mendonga transformation, is described 
by the dual Rankine-Hugoniot conditions (62-64). 

Since the coordinates x and t are continuous on S, and crucially because the 
relations between the currents are linear, the conditions (62-64) are equivalent to 
the jump conditions obtained from mass, momentum and energy conservation, 
in (42-44). In fact, these two sets of equations imply, and are implied by, each 
other. In conclusion, the dual set of jump conditions associated with mass, 
boosts and expansions, is completely equivalent to the usual Rankine-Hugoniot 
conditions and may be used for an independent description of the shock. 

6 The Entropy Condition 

For a polytropic gas, by choosing e = xp 70 , we can rewrite Eq. (5) as 

D X = (65) 

In we defined an isentropic flow to be one for which \ — constant. For a 
general flow, it followed that \ transforms like a scalar. For a polytropic gas, it 
is also well-known [H] that x is related to the specific entropy (entropy per unit 
mass), S as follows: 

S - So = C v \og[x{ P Vy°] (66) 

where C v = R/("fo ~ l)j R being the universal gas constant divided by the 
molecular weight of the particular gas, V the volume and So an appropriate 
constant. It is obvious from this equation that as a particle of the medium moves 
about, the specific entropy at the moving particle remains constant under an 
SL(2, R) transformation. Hence, under an SL(2, R) transformation, a physical 
shock gets mapped to a physical shock. 

We now require the positivity of viscosity to be preserved under an SL(2, R) 
transformation - a requirement that guarantees that the shock respects the 
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entropy condition. As already pointed out (see eq. (12)), in three-dimensional 
space, the viscosity fields transform as follows: 

rf = (ft + sfri and C'^^t + dfc (67) 

Thus the transformation properties of the viscosity fields are similar to p i.e., 
they transform like scalar densities. Hence, if we integrate the viscosity field over 
all space, to get the viscosity, it is an invariant under the SX(2, R) transforma- 
tions. Likewise, the specific viscosity (viscosity per unit mass), is an invariant. It 
follows that the positivity of the viscosity is maintained without any additional 
restrictions on the SL(2, R) parameters. 

7 Conclusions 

In this paper, we extended the analysis of jo] to discuss the symmetries of 
discontinuous flows in fluid dynamics. The maximal kinematical invariance 
group of an ideal, polytropic fluid is Q = SL(2, R) A G, not just for smooth, but 
for discontinuous flows also. This is made manifest by writing the fluid equations 
in their conservation law form. New conservation laws follow from a direct 
application of Noether's theorem, enabling us to construct a dual set of Rankine- 
Hugoniot shock conditions. The SL(2, R) transformations map the standard 
Rankinc-Hugoniot shock conditions to the dual ones and vice versa. These 
transformations also respect the entropy conditions that physical shocks need 
to satisfy. Hence we conclude that, under these transformations, an explosion 
gets mapped to an implosion, thus offering a theoretical explanation for the 
intriguing observations of ^|2IH]- 
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